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ABSTRACT 

We calculate the Isgur-Wise function by measuring the elastic scattering amplitude of a D meson in the quenched 
approximation on a 24 3 x 48 lattice at (3 = 6.2, using an 0(a)-improved fermion action. We use this result, in 
conjunction with heavy-quark symmetry, to extract \V c b\ from the experimentally measured B — > D*li> differential 
decay width. 



CO 



Heavy-quark symmetry enables all the non-pertur- 
bative, strong-interaction physics for semi-leptonic B — > 
D, D* decays to be parametrized in terms of a sin- 
gle universal function, £(u>), of u> = v ■ v', where v 
and v' are the four- velocities of the B and D mesons 
respectively [1, 2]. is known as the Isgur-Wise func- 



tion and is normalized at zero recoil: £(1) = 1[2]. With 
£(w), one can extract the Cabibbo-Kobayashi-Maskawa 



^matrix element V c b from experimental measurements of 
/■vthe rate for B — > D*li> decays [3]. We report here on a 
"lattice QCD calculation of the Isgur-Wise function and 
OOon the corresponding determination of V c b- We also 
perform a qualitative test of the flavor component of 
heavy-quark symmetry. 

CNI To obtain the Isgur-Wise function, we evaluate the 
<Nelast 

ic scattering matrix element (D(p')\cj fl c\D(p)) on 
^-^the mass shell [4]. Because the electromagnetic current 
Q\cj^c is conserved, this matrix element can be parame- 
"^ trized in terms of a single form factor: 

Qh (D(p')\c^c\D(p)) =m D (v + v'yh el {Lo) , (1) 

ftidiere p^ = mov^ and to = v ■ v' is the four- velocity 
irecoil. In the limit of exact heavy-quark symmetry this 
form factor is simply £(u>). 

There are two sources of corrections to this simple 
result: 



h el (w)= [l + /3 c 'M+ 7 e 'M]£M 



(2) 



The first correction, (3 el (ui), stems from perturbative 
QCD corrections to the heavy-quark current. We ob- 
tain this correction from Neubert's short distance ex- 
pansion of heavy-quark currents [5]. The second correc- 
tion, 7 e '(w), is due to higher-dimension operators with 
coefficients proportional to inverse powers of the charm 
quark mass. 7 e '(o;) is difficult to quantify because it in- 
volves the light degrees of freedom and is therefore non- 
perturbative. Luke's theorem[6], however, guarantees 
that there is no (AQCD/(2m c )) correction to h el (u>) 
at zero recoil. Moreover, model estimates of this cor- 
rection appears to remain well below 3% over the range 
of experimentally accessible recoils [7]. That this correc- 
tion is small is corroborated by the fact that we see no 



differences in the ratio h el (u>) / (1 + /3 el (ui)) computed for 
two values of the heavy-quark mass. Thus, we will ne- 
glect 7 e '(w) in extracting the Isgur-Wise function from 
h el (u). As defined in Eq. (2), £(u>) is renormalization- 
group invariant and normalized to one at u> = 1[5]. 

We work in the quenched approximation on a 24 3 x 
48 lattice at (3 = 6.2, which corresponds to an inverse 
lattice spacing a -1 = 2.73(5) GeV, as determined from 
the string tension [8]. Our calculation is performed on 
sixty SU (3) gauge field configurations (for details see 
Ref. [8]). The mesons are composed of a propagat- 
ing heavy quark with a mass around that of the charm 
quark, and a light antiquark with a mass around that 
of the strange quark. To reduce discretization errors, 
the quark propagators are calculated using an 0(a)- 
improved act ion [9]. This improvement is particularly 
important here since we are studying the propagation 
of quarks whose bare masses are around one half the 
inverse lattice spacing. Our statistical errors are calcu- 
lated according to the bootstrap procedure described in 
Ref. [8]. 

To obtain the matrix element (D(p')|c7 M c|D(p)), we 
calculate the ratio of three-point correlators, 



= Ex, y e-^+»'^(J(y)V»(x) Jt(0)) 
[X> ~ Z x , y (J(v)V°(x)Jl(0)) 



(3) 



where J is a spatially-extended interpolating field for 
the D meson[10], V M is the 0(a)-improved version of 
the vector current C7 M c[ll] and p=p' + q. To evaluate 
these correlators, we use the standard source method[12]. 
We choose t y = 24 and symmetrize the correlators 
about that point using Euclidean time reversal. We 
evaluate for three values of the light-quark mass 
(ki = 0.14144, 0.14226, 0.14262) which straddle the 
strange quark mass (given by k s = 0.1419(1)[13]); two 
values of the heavy-quark mass (kh = 0.121,0.129) 
around that of the charm quark (given by re c ~0.129[14]); 
and several values of the initial and final .D-meson mo- 
menta, all less than {it /\2a)\/2. 
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Figure 1: /i e '(w)/(l + /3 e '(w)) vs w. 



Figure 2: vs cj. 



In the limit t x , t v — t x — > oo, we fit to 



m D Z(p 2 )Z(p' 2 ) 
2EE' Z 2 (0) 



(S-S'Jt^ + CS'-mr))^ 



<D(p')|c7°c(0)|D(p)> , (4) 



where E (E') is the energy of the initial (final) D me- 
son and Z(p 2 ) = (0| J(0)\D(p)). We see a plateau in 
R about t x = 12, and fit R to a constant for t x = 
11, 12, 13. Multiplying this constant by suitable wave- 
function and energy factors, obtained from fits to two- 
point functions, we get (D(p')\cj°c\D(p)) . 

The ratio h el (co)/(l + (3 el (co)) that we obtain for 
K[ = 0.14144, the heaviest of our light quarks, is plot- 
ted in Fig. 1. Data for both values of the heavy-quark 
mass, m/j, are shown. The squares correspond to nihi = 
1.4 GeV (kh = 0.129) and the octogons to rrih2 = 1-9 
GeV (kh = 0.121)*. The two sets of data agree very 
well. It is difficult, however, to use this information to 
quantify the 1 jm c correction since our range of heavy- 
quark masses is small. Only (m^ — mhi)/mh2 — 1/4 of 
the full l/m c correction can be resolved. The fact that 
we see no sign of this correction indicates that it must 
be smaller than four times a few percent. In Fig. 1, the 
solid curve is a two-parameter fit to s£ p (u>), where £, p (ui) 
is Stech's relativistic-oscillator parametrization[7]: 



*>) = ^exp (-(2^-1)^) 



(5) 



and p 2 = —£' (1). The parameter s was added to absorb 
uncertainties in the overall normalization of our data. 



* m /il(/i2) i s the spin averaged mass of pseudoscalar and vec- 
tor mesons composed of a heavy-quark, /ii(/i2), and a massless 
antiquark[14], less 0.5 GeV (p. 79 of Ref. [15]) 



(T = 1.5 ±2 and s = 0.95 +J with a x 2 / do / = 
0.8. Other parametrizations for give very similar 
results. 

If we fit our data to instead of s£ p (cu) (an 

equally valid procedure, in principle, for determining 
p 2 ), we find p 2 = 2.1 with a x 2 /Aof = 2.9 (dotted 
curve in Fig. 1). To accommodate the spread in values 
for p 2 given by our two procedures, we assign errors to 
p 2 which encompass all values consistent with both pro- 
cedures. These errors include systematic uncertainties, 
but only to the extent that the deviation of s from 1 , in 
our first fit, is an indication of systematic errors. The 
central value we choose for p 2 is the one given by our 
first fit since this fit is designed to absorb uncertain- 
ties in the overall normalization of our data. Thus, for 
0.14144 we quote p 2 



Ki 



In Fig. 2 we plot the results obtained from a covari- 
ant and linear extrapolation of our data for h el (uj)/(l + 
/? e '(w)) for the three values of the light-quark mass to 
Ki = Kent- These results correspond to a meson com- 
posed of a charm quark (kh = 0.129) and a massless 
antiquark. They are our Isgur-Wise function. With a 
two parameter fit to s£ p (lo), we get p 2 = 1.2 lg and s = 

0.94 ±1 with a xV do / = °- 9 ( solid line in Fi S- 2 ). Forc- 
ing s to be 1, we find p 2 = 1.7 1% with a x 2 /do/ = 1.6 
(dotted line in Fig. 2). Using the same procedure as the 
one used to determine the errors and central value for 
p 2 when Ki = 0.14144, we obtain, as our best estimate 
for p 2 , when K; = K cr i t : 



1.2 



+ 7 



(6) 



Dependence of on light-quark mass is mild. Fur- 
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Figure 3: |V c &|£(k>) vs to. The solid line is |V c b|£p(w) 
with p 2 = 1.2. 



thermore, our result for p 2 agrees with most other de- 
terminations of this parameter[3, 16], apart from the 
sum-rule result of Ref. [17] which lies below our error 
bars. In particular, our result for p 2 agrees with the 
very recent lattice result of Ref. [18] obtained with Wil- 
son fermions, although the details and systematics of 
the two calculations are different. 

We can now obtain |V C &|. Heavy-quark symmetry 
enables one to write the differential decay rate for B — > 
D*lv decays in terms of the Isgur-Wise function: 



dT 

duo 



« \v cb \ 2 a^) 2 , 



(7) 



where the proportionality factors are described in Ref. [3] . 
Thus, experimental data for dT/dco yield a measure- 
ment of |V c &|£(k>) versus lo. To get |V C &|, we fit |V c b|£p(w) 
to these data with £ p (u>) given by Eq. (5) and p 2 fixed 
by our lattice calculation to the value given in Eq. (6). 
In Fig. 3 we show this least-x 2 -fit for new ARGUS 
data[19]. The resulting value for |V c f,| is 



IK 



48ps 



0.038 



+2+8 



(8) 



with a x 2 /do/ = 1.1. The same fit to the weighted 
average of older CLEO and ARGUS data[20] gives |V c6 | 

^/r B /1.48ps = 0.036 1 2 .^ with a xV do / = °- 6 - In 
both cases, the first set of errors is due to experimental 
uncertainties, while the second is due to the uncertainty 
in our lattice determination of p 2 . The B-meson lifetime 
used above is the central value of the lifetime, t b = 
1.48(10)ps, quoted in Ref. [21]. 
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